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Abstract. We derive and analyze monotone difference-quadrature schemes 
for Bellman equations of controlled Levy (jump-diffusion) processes. These 
equations are fully non-linear, degenerate parabolic integro-PDEs interpreted 
in the sense of viscosity solutions. We propose new "direct" discretizations of 
the non-local part of the equation that give rise to monotone schemes capable 
of handling singular Levy measures. Furthermore, we develop a new general 
theory for deriving error estimates for approximate solutions of integro-PDEs, 
which thereafter is applied to the proposed difference-quadrature schemes. 
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1. Introduction 

In this article we derive and analyze numerical schemes for fully non-linear, 
degenerate parabolic integro partial differential equations (IPDEs) of Bellman type. 
To be precise, we consider the initial value problem 



+ sup { - L a [u]{t,x) + c a {t,x)u~ f a (t,x) - J a [u]{t lX )\ = in Q T , (1.1) 

u(0,x)=g(x) inR N , (1.2) 
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where Q T := (0, T] x R N and 

L Q [0](i,x) := tr[a a {t,x)D 2 (j)] + b a (t,x)D(j), 

J a [<j>](t,x) := / (0(t,a; + » ? a (i,a;,z))-0-l|,|<it ? Q (t J a:,*)D0)i/(<te), 
Jr m \{o} v 7 

for smooth bounded functions 0. Equation (jl.lj) is convex and non-local. The 
coefficients a a ,i] a ,b a ,c a , f a ,g are given functions taking values respectively in S w 
(N x N symmetric matrices), R N , R N , M, K, and K. The Levy measure i/(dz) is 
a positive, possibly singular, Radon measure on R M \{0}; precise assumptions will 
be given later. 

The non-local operators J a can be pseudo-differential operators. Specifying 
rj = z and v{dz) = \ z \n+-, dz, 7 £ (0, 2), give rise to the fractional Laplace operator 

J = (—A) 7 / 2 . These operators are allowed to degenerate since we allow 77 = 
for z ^ 0. The second order differential operator L a is also allowed to degenerate 
since we only assume that the diffusion matrix a a is nonnegative definite. Due to 
these two types of degeneracies, equation (II. ip is degenerate parabolic and there is no 
(global) smoothing of solutions in this problem (neither "Laplacian" nor "fractional 
Laplacian" smoothing) . Therefore equation (jl.lj) will have no classical solutions in 
general. From the type non-linearity and degeneracy present in (jl.lj) the natural 
type of weak solutions are the viscosity solutions [SO] [25] . For a precise definition 
of viscosity solution of (jl.ip we refer to [57]. In this paper we will work with 
Holder/Lipschitz continuous viscosity solution of (|1.1[) - (|1.2[) . For other works on 
viscosity solutions and IPDEs of second order, we refer to [21 HI [3 HI EH E3 H3 
[23 E7J HO] and references therein. 

Nonlocal equations such as (II. 1|) appear as the dynamic programming equation 
associated with optimal control of jump-diffusion processes over a finite time horizon 
(see [37J [35] [T5]). Examples of such control problems include various portfolio 
optimization problems in mathematical finance where the risky assets are driven 
by Levy processes. The linear pricing equations for European and Asian options in 
Levy markets are also of the form (|l.ip if we take A to be a singleton. For more 
information on pricing theory and its relation to IPDEs we refer to [18] . 

For most nonlinear problems like (jl.ip ~ p.2p . solutions must be computed by a 
numerical scheme. The construction and analysis of numerical schemes for nonlinear 
IPDEs is a relatively new area of research. Compared to the PDE case, there are 
currently only a few works available. Moreover, it is difficult to prove that such 
schemes converge to the correct (viscosity) solution. In the literature there are two 
main strategies for the discretization the non-local term in (jl.ip . One is indirect 
in the sense that the Levy measure is first truncated to obtain a finite measure 
and then the corresponding finite integral term is approximated by a quadrature 
rule. Regarding this strategy, we refer to [THl HH] (linear or obstacle problems) 
and (29] [16] (general non- linear problems). The other approach is to discretize the 
integral term directly. Now there are 3 different cases to consider depending on 
whether (i) S\ z \<i u {dz) < 00, (ii) S\ z \<i \ z W(dz) < 00, or (iii) S\z\<i \ z \ 2 v(dz) < 00. 
Case (i) is the simplest one and has been considered by many authors, see, e.g., 
[33 QI1 El [55] [3 [55] and references therein. Case (ii) was considered in [U [33 [53], 
and case (iii) in [33 [53]. Most of the cited papers restrict their attention to linear, 
non-degenerate, one-dimensional equations or obstacle problems for such equations. 

One of the contributions of this paper is a class of direct approximations of the 
non-local part of (jl.ip , giving rise to new monotone schemes that are capable of han- 
dling singular Levy measures and moreover are supported by a theoretical analysis. 
The proposed schemes are new also in the linear case. As in [T] (cf. also [33] for a 
related approach), the underlying idea is to perform integration by parts to obtain a 
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bounded "Levy" measure and an integrand involving derivatives of the unknown so- 
lution. In pP, one-dimensional, constant coefficients, linear equations (and obstacle 
problems) are discretized under the assumption j\ z \ <x \z\ v(dz) < oo. Their schemes 
are high-order and non-monotone, but not supported by rigorous stability and con- 
vergence results. In this paper we discretize general non-linear, multi-dimensional, 
non-local equations without any additional restrictive intcgrability condition on the 
Levy measure. More precisely, we provide monotone difference-quadrature schemes 
for (II. I j) - (|l ,2p and prove under weak assumptions that these schemes converge with 
a rate to the exact viscosity solution of the underlying IPDE. The schemes we put 
forward and our convergence results apply in much more general situations than 
those previously treated in the literature. 

The second main contribution of this paper is a theory of error estimates for a 
class of monotone approximations schemes for the initial value problem (jl.f [) - (|1.2p . 
We use this theory to derive error estimates for the proposed numerical schemes. For 
IPDEs in general and non-linear IPDEs in particular, there are few error estimates 
available, see [T9l [55] for linear equations and [221 [HI US] for non-linear equations. 

Error estimates involving viscosity solutions first appeared in 1984 for first order 
PDEs [UJ, in 1997/2000 for convex 2nd order PDEs [13 HQ, and in 2005/2008 for 
IPDEs [T§1 123] • The results obtained for IPDEs, including those in this paper, are 
extensions of the results known for convex second order PDEs, which are based on 
Krylov's method of shaking the coefficients [31] . Krylov's method produces smooth 
approximate subsolutions of the equation (or scheme) that, via classical comparison 
and consistency arguments, imply one-sided error estimates. Based on this idea, 
there are currently two types of error estimates for convex second order PDEs: (i) 
optimal rates applying to specific schemes and equations (cf., e.g., [H21 [2S]) and 
(ii) sub-optimal rates that apply to "any" monotone consistent approximation (cf., 
e.g., [51[3T]). In particular, type (i) results apply when you have a priori regularity 
results for the scheme, while type (ii) results do not require this. 

In this paper we provide error estimates of type (ii) , whereas earlier results for 
IPDEs are of type (i), see [EH [23 HH EI] ■ The problem with type (i) results is the 
difficulty in establishing the required priori regularity estimates. In the PDE case 
this can be achieved for particular schemes [32] [26 , and attempts to generalize these 
schemes to the IPDE setting have only been partially successful [TTJ [H] , since the 
required regularity estimates have been obtained only through unnaturally strong 
restrictions on the non-local terms. In [16] the Levy measure is bounded and in 
[iT] the Levy measure is either bounded or the integral term is independent of x 
with an (essentially) one-dimensional Levy measure. Of course, by a truncation 
procedure only bounded Levy measures need to be considered [TBI US] , but such 
approximations may not be accurate and the resulting error estimates blow up as 
the truncation parameter tends to zero. An advantage of the error estimates in the 
present paper is that they apply without any such restrictions. In particular, we 
can handle naturally any singular Levy measures directly in our framework. 

To prove our results we extend the approach of [5] to the non-local setting. To 
this end, we have to invoke a switching system approximation of (|1.1|) (see Section 
[SJ). Switching systems of this generality have not been studied before. In paper 
[12j . we provide well-posedness, regularity, and continuous dependence results for 
such systems. We also prove that the value function of a combined switching and 
continuous control problem solve the switching system under consideration. 

The remaining part of this paper is organized as follows: First of all, we shall 
end this introduction by listing some relevant notation. In Section [2] we list a 
few standing assumptions and provide corresponding well-posedness and regularity 
results for the IPDE problem (jl . 1(1 - (|1 .2[) . In Section [3] we present a rather general 
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approximation scheme for this problem, and show that it is consistent, monotone, 
and convergent. Error estimates for general monotone approximation schemes are 
stated in Section |4] In Section [5] we present new direct discretizations of the non- 
local term in and prove that these discretizations are consistent, monotone, 
and also satisfy the requirements introduced in Section [3] The switching system 
approximation of (|l.ip is introduced and analyzed in Section [HI The obtained 
results are utilized in Section[7]to prove the error estimate stated SectionHl Finally, 
in Appendix E] we give a standard example of a (monotone) discretization of the 
local PDE part of (jl.ip that satisfies the requirements of Section [3] 

We now introduce the notation we will use in this paper. By C, K we mean 
various constants which may change from line to line. The Euclidean norm on any 
R d -type space is denoted by | • |. For any subset QcRx M. N and for any bounded, 
possibly vector valued, function on Q, we define the following norms, 

... \w(t,x) — w(t,y)\ 
\w\o-= sup \w(t,x)\, Mi = Mo+ SU P i a — ; 

(t,x)eQ (t,x)^(s, y ) \t-s\i +\x-y\ 

Note that if w is independent of t, then is the Lipschitz (or W 1,00 ) norm of w. 
We use Cb{Q) to denote the space of bounded continuous real valued functions on 
Q. Let p(t,x) be a smooth and non-negative function on R x with unit mass 
and support in {0 < t < 1} x {\x\ < 1}. For any e > 0, we define the mollifier p e by 

Pe{t,x) := -jjTnp(-^,-)- (1-3) 



In this paper we denote by h the vector 

h = {At, Ax, Az) > 0, 

and any dependence on At, Ax, or Az will be denoted by subscript h. The grid is 
denoted by Qh and is a subset of Qt which need not be uniform or even discrete in 
general. We also set Q% = Q h n {t = 0} and Q+ = Q h n {t > 0}. 

2. WELL-POSEDNESS & REGULARITY RESULTS FOR THE BELLMAN EQUATION 

In this section we give some relevant well-posedness and regularity results for the 
Bellman equation (|l.l|) - l)1.2)l . To this end, we impose the following assumptions: 

(A.l) The control set A is a separable metric space. For any a £ A, a a = 
^a a a aT , and cr a ,b a , c a , / Q , rj a are continuous in a for all x, t, z. 

(A. 2) There is a positive constant K such that for all a. € A, 

Mi + Kii + mi + Hi + irii <k. 

(A. 3) For every a € A and z € R M there is an A > such that 

le-^VO,-^)!! <K(\z\M) and |e" A IV(t, x , -)|i < K. 
(A. 4) v is a positive Radon measure on R M \ {0} satisfying 

f \z\ 2 v{dz)+ f e (A+e)|z| Kdz) < K 



for some K > 0, e > where A is defined in (A. 3) 



Sometimes we need the following stronger assumptions than (A. 3) and (A. 4) 
(A. 4') v is a positive Radon measure having a density k{z) satisfying 

0<fc(^)< e | z ( |M+T' foral1 *eR M \{0}, 



for some 7 <E (0, 2), e > 0, where A is defined in (A. 3) 
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(A. 5) Assume that (A. 3) holds and let 7 as in (A. 4'). There is a constant K such 
that for every a £ A and z £ M M 



for all 



|£>^ a (v^)|o + |^( v ^)| <Xe A W, 

k = I = 1 when 7 = 0, 

Me {1,2} when 7 £ (0,1), 

k £ {1, 2, 3, 4}, J € {1, 2} when 7 € [1, 2). 



Assumptions (A.l) (A.4) are standard and general. The assumptions on the 
non-local term are motivated by applications in finance. Almost all Levy models 
in finance are covered by these assumptions. It is easy to modify the results in 
this paper so that they apply to IPDEs under different assumptions on the Levy 
measures, e.g., to IPDEs of fractional Laplace type where there is no exponential 
decay of the Levy measure at infinity. Finally, assumption (A. 5) is not strictly 
speaking needed in this paper. We use it in some results because it simplifies some 
of our error estimates. 

Under these assumptions the following results hold: 



Proposition 2.1. Assume (A.l) - (A.4) 



(a) There exists a unique bounded viscosity solution u of the initial value problem 
(jl.ip - (|1.2p satisfying \u\\ < 00. 

(b) If u\ andu2 are respectively viscosity sub and supersolutions of (jl.ip satisfying 
«i(0, •) < ^2(0, •); then u\ <U2- 

The precise definition of viscosity solutions for the non-local problem (|l.l|) - ljl.2|) 
and the proof of Proposition 1 2 . 1 1 can be found in [37j, for example. 

3. Difference-Quadrature schemes for the Bellman equation 

Now we explain how to discretize (|l.ip - (|1.2[) by convergent monotone schemes 
on a uniform grid (for simplicity). We start by the spatial part and approximate 
the non-local part J a as explained later in Section \E\ and the local PDE part L a by 
a standard monotone scheme (cf. [33] and Appendix EJ ■ The result is a system of 
ODEs in AxZ N x (0,T): 

u t + sup { - L%[u](t,x) + c a (t,x)u- f a {t,x) - J%[u](t,x)\ = 0, 

where Lh and Jh are monotone, consistent approximations of L and J, respectively. 

Then we discretize in the time variable using two separate ^-methods, one for 
the differential part and one for the integral part. For ?9, £ [0, 1], the fully discrete 
scheme reads 



U n = U-- 1 - At sup { - OLlPTp - (1 - 6)Ll[U]f x + c^U^ 1 (3.1) 

_ f ^-X_ §J a [u] n_ il _^ja [u] n-l^ j n g. 



Uf 3 =g{xp) in Q° h , (3.2) 

where Q h = AxZ N x At{0, 1, 2, . . . , £} and U% = U(t n ,xp), = f a (t n ,x p ), 
etc., for t n = nAt (n £ N ) and xp = [3Ax {/3 £ 1 N ). 

The approximations Lh and Jh are consistent, satisfying 

|L Q H - L?M\ < K L (\D 2 <f>\ Ax + l^loAx 2 ), (3.3) 

\{\D 2 (j)\o + \D 4 (f)\o) when 7= [1,2), 
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for smooth bounded functions <f> an d where 7 £ (0, 2) is defined in (A. 4'). They are 
also monotone in the sense that they can be written as 



L%[4>](t n ,Xp) = W*"' X P) ~ 0(*~> x fi)] 



ith I 



>0, (3.5) 



Jh\<t>]{tn>Xp) = 53 Cj^*"'^' ~<l>(tn,Xp)] with j°^>0, (3.6) 

for any (3 £ Z N and n £ No- We also assume without loss of generality that 
J-'fl = ® = t'p B f° r au ^ e ^ N ■ ^ ne sum P-^P is always finite, while the sum 
(|3 6p is finite if the Levy measure v is compactly supported. With 7 S [0, 2) defined 
in (A. 4') and Ax < 1, we also have that 



f /3 



•a,n 



& - ^ sup 



pel 



upf |a Q | Aa;" 



\b a \ Ax 



E .o,n <i( - a -1 



(3.7) 
(3.8) 



From (|3.3p and (|3.4[) it immediately follows that the scheme (|3.1[) is a consistent 
approximation of fTTTj) , with the truncation error bounded by 



i|0 tt | o A< + sup ||L Q [C - L«H"| + |J a [0] n - J^CIo 

+ (1 - 9)\L a [<^- x - L Q [0]"| O + (1 - 0)| rW- 1 - J Q [0]"| O }, 
for smooth functions 4>. The last two terms are again bounded by 

At sup {|L a [&]| + k Q [0t]|o} < KAt[\d t Dcj>\ + \d t D 2 cj>\ o y 



(3.9) 



(3.10) 



Under a CFL condition, the scheme 



is also monotone, meaning that there 



are numbers b ~ (a) > such that it can be written as 



for all (xg,t n ) £ Q]\. From ()3.5|) and ()3.6p . we see that 



b P-,P («) = 



1 + At6lf m + AW~j a ' m 
I- At 



a,m a,m 



in,m/ \ 
b 0,3 (") 



a.m 



when m = n, 
when 771 = 71 — 1 , 

when 777 = 77, 
when 777 = 77 — 1 , 



where (3 ^ (3 and other choices of m give zero. These coefficients are positive 
provided the following CFL condition holds: 



At 



(1-0)1/9 ~ C /3 



< 1 for all a, f3, m, 



(3.12) 



or alternatively by (|32J) and ([3"%]) . if a Q ^ 0, c" > 0, Ax < 1, 



At 



(1 - 6)KiCAx~ 2 + (1 - &)KjAx- 



< 1. 



Existence, uniqueness, and convergence results for the above approximation 
scheme are collected in the next theorem, while error estimates are postponed to 
Theorem [O] in Section g) 
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Theorem 3.1. At 



JA~T]\\(AI]\ (A4'), 



(GEU), and igUl) . 



(cij There exists a unique bounded solution Uh of (|3.1[) - (|3.2p . 

(7>j The scheme is L°° -stable, i.e. \U h \ < e sup ° l c °l *" [\g\ + t n sup Q |/ Q | ] . 

(c) Uh converge uniformly to the viscosity solution u of (|l.lj) - (|1.2|) as ft, — > 0. 

Proof. The existence and uniqueness of bounded solutions follow by an induction 
argument. Consider t — t n and assume C/' i_1 is a given bounded function. For 
e > we define the operator T : U n — > J7" by 



rjra 



e • (left hand side of ([5TTT]) ) for all /? e Z 



Note that the fixed point equation {/" = TU" is equivalent to equation 
Moreover, for sufficiently small e, T is a contraction operator on the Banach space 
of bounded functions on Ax Z N under the sup-norm. Existence and uniqueness 
then follows from the fixed point theorem (for U n ) and for all of U by induction 
since U° = g\go is bounded. 

To see that T is a contraction we use the definition and sign of the ^-coefficients: 



< sup • 



{ [l - e[l + At(8l a p n + (L7£ - U$) + eAt(9l a / 1 + tij™' n )\U n - U n \ } 



< {l-e)\U % 



U r ' 



provided 1 — s(l + At(9l^' n + > for all a,(3,n. Taking the supremum 

over all (3 and interchanging the role of U and U proves that T is a contraction. 

Much the same argument, utilizing (|3.11|) . establishes that Uh is bounded by a 
constant independent of h: 



|f/"|o < (l + Aisup|c Q | r |5|o + nAtsup|r|o 



< gSUPo |c"|o*r> 



|.g|o + i„sup|/ a | 



In view of this bound, the convergence of Uh to the solution u of (ll.ll) - (|i.2p follows 
by adapting the Barles-Souganidis argument [S] to the present non-local context. 
Alternatively, convergence follows from Theorem 14.31 if we also assume (A. 5). □ 



Remark 3.1. 

a. One suitable choice of J% will be derived in Section El while for there 
are several choices that satisfies (|3.3p and Q3.5p . e.g., the scheme by Bonnans and 
Zidani [13 or the (standard) schemes of Kushner [33]. In Appendix lAlwe show that 
one of the schemes of Kushner fall into our framework if a a is diagonally dominant. 

b. For the differential part, the choices 9 — 0, 1, and 1/2 give explicit, implicit, 
and Crank-Nicholson discretizations. When -d > 0, the integral term is evaluated 
implicitly. This leads to linear systems with full matrices and is not used much in 
the literature. 

c. By parabolic regularity D 2 ~ <9 t and (|3.10p is similar to At\4>tt\a- When 
9 = 1/2 = i9 the scheme (I3.ip (Crank-Nicholson!) is second order in time 0(At 2 ) 
and (|3.9p is no longer optimal. 

d. When 7 = the leading error term in Jh[u] (see (|3.4p ) comes from difference 
approximation of the term Du J r\v. This difference approximation also give rise 
the term Aa;" 1 in ([3T8]) . 
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4. Error estimates for general monotone approximations 

In this section we present error estimates for nonlinear general monotone ap- 
proximation schemes for IPDEs. As a corollary we obtain an error estimate for the 
scheme (|3.1I) - (|3.2I) defined in Section [31 These results, which extend those in [5] to 
the non-local IPDE context, can be applied to "any" Levy-type integro operator. 
Earlier results apply to either linear problems, specific schemes, or restricted types 
of Levy operators, see (33 EH US HI]- In particular, previous error estimates do 
not apply to the approximation scheme (|3.ip . 

Let us write as u t + F[u] = where F[u] := F(t,x,u,Du,D 2 u,u(t,-)) 

denotes the sup part of (|l.ip . We write approximations of u t + F[u] — as 

S{h,t,x,u h (t,x), [u h ]t,x) = in (4.1) 
u h (0,x)=g h (x) in G° h , (4.2) 

where S is the approximation of defined on the mesh Qh C Qt with "mesh" 
parameter h = (At, Ax, Az) (time, space, quadrature parameters). The solution is 
typified by Uh and by [tijt.j; we denote a function defined at (t,x) in terms of the 
values taken by ut evaluated at points other than (t,x). Note that the grid does 
not have to be uniform or even discrete. 

We assume that (|4.1j) satisfies the following set of (very weak) assumptions: 

(51) (Monotonicity) There exist A, /i > 0, ho > such that, if \h\ < ho, u < v 
are functions in Cb(Qh) and 4>(t) = e M *(a + bt) + c for a,b,c> 0, then 

S(h,t,x,r + 4>(t), [u + <t>] tjX ) > S(h,t,x,r, [v] t , x ) + - - Ac in 0+. 

(52) (Regularity) For each h and <fr G Cb(Qh), the mapping 

(t,x) h-> S(h,t,x,4>(t,x), [4>}t,x] 

is bounded and continuous in and the function r i— > S(h, t, x, r, \<$\t.x) is 
uniformly continuous for bounded r, uniformly in t, x. 

(53) (i) (Sub-consistency) There exists a function E\(K ,h,e) such that, for 
any sequence {</> e } e of smooth bounded functions satisfying 

\df° D p ' 4> t \ < # e i- 2 AH/3'l in q Tj for any ^ e N,/3' e N w , 
where |/3'| = YliLi P'n the following inequality holds in Q~£: 
S(h, t, x, e (t, x), [<t> e ] t , x ) < c>>et + F(t, x, ^,D^,D 2 <j>, <p e (t, •)) + Ei{K, h, e). 

(S3) (ii) (Super-consistency) There exists a function E2(K,h,e) such that, 
for any sequence {4>e}e of smooth bounded functions satisfying 

\d^° D f3 ' (t) e \ < if e l-2/3o-|/3'l m Q T; for any p Q G N)/ g/ g 

the following inequality holds in Q^: 

S(h, t, x, <f> e (t, x), [0 e ] M ) > d f A + F(t, x, ^,D<p t ,D 2 <j), <j> e (t, •)) - E 2 (K, h, e). 

Remark 4.1. In (S3), we typically take e = w e * p e for some sequence (w e ) e 
of uniformly bounded and Lipschitz continuous functions, and p e is the mollificr 
defined in Section [T] 

Remark 4.2. Assumption (SI) implies monotonicity in [u] (take <fi — 0), and 
parabolicity of the scheme (|4.ip (take u — v). This last point is easier to understand 
from the following more restrictive assumption: 
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(SI') (Monotonicity) There exist A > 0, K > such that if u < v;u, v £ C b (Qh) 
and 4> : [0, T] -*1 smooth, then 

S(h, t,x,r + <p(t), [u + (f>] tiX ) 

> S(h,t,x,r, [v] t<x ) + <f>'(t) - KAt\4>"{t)\ - A0+(t). 

It is easy to see that (SI') implies (SI), cf. [8]. 

The main consequence of (SI) and (S2) is the following comparison principle 
satisfied by scheme (|4.1I) (for a proof, cf. [5]): 

Lemma 4.1. Assume (SI), (S2), 51,32 G Cb{Gh), andu,v e Cb{Qh) satisfy 

S{h,t,x,u(t,x),[u]t, x ) < gi and S(h, t, x, v(t, x), [v] t , x ) > gi in . 
Then, for A and \x as in (SI), 

u-v< e"*|(u(0, •) - w(0, -)) + |o + 2te t |(«? 1 - g 2 )+\ . 
The following theorem is our first main result. 



Theorem 4.2 (Error Estimate). Assume (A.l)- (A4), (SI), (S2) hold, and that 
the approximation scheme (14.1|) - (|4.2p has a unique solution Uh G Cf,(Gh)> for each 
sufficiently small h. Let u be the exact solution of (|l.lj) - (|1.2| . 

a) (Upper Bound) If (S3)(i) holds, then there exists a constant C, depending 
only on IJ,,K in (SI) and (A. 2) such that 

u-u h < e^Hg-, g h ) + \o + Cmin (e + /i,e)) in Q h . 

b) (Lower Bound) If (S3)(ii) holds, then there exists a constant C , depending 
only on IJ,,K in (SI) and (A. 2), such that 

U - u h > -e^Kg - g h y\{) - Cmin (e^ + E 2 {\u\\, h, e)) in Q h . 

£>0 

We prove this theorem in Section [7J 



Remark 4.3. Theorem 14.21 applies to all Levy type non-local operators. Note 
that the lower bound is worse than the upper bound, and may not be optimal. In 
certain special cases it is possible to prove better bounds, however until now such 
results could only be obtained in the non-degenerate linear case [351 UH] or under 
very strong restrictions on the non-local term [111 116] . More information on such 
non-symmetric error bounds can be found in [5]. 

Remark 4.4. For a finite difference-quadrature type discretization of (jl.lj) . the 
truncation error would typically look like 

\<t> t + F(t, x, 0, D(j>, D 2 <j), <t>(t, •)) - S(h, t, x, 0(t, x), [<t>] t , x )\ 

< KJ2 \df°D f3 '°<j>\ At k v« +Kj2\df 1 D (3 'i ( f)\( ) Ax k ^ + Kj^ \df 2 D^(f>\ Az k ^ , 

00 01 02 

where f3 = (Po,f3' ), ft. = ft = (P 2 ,(3' 2 ) are multi-indices and k0 o ,kp 2 ,kp 2 

are real numbers. In this case, the function E in (S3) is obtained by taking <f> := <fi e 
in the above inequality: 

E 1 =E 2 = KK [e l - 2<3 "- ]0 '°\ At k ?° + e 1 - 2 ^"^! Ax k ^ + e^^-^AzM . 

00,01,02 

An optimization with respect to s yields the final convergence rate. Observe that 
the obtained rate reflects a potential lack of smoothness in the solution. 

We shall now use Theorem 14.21 to prove error estimates for the finite difference- 
quadrature scheme 
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Theorem 4.3. Assume\(AlJW(A. 3)\ (A. 4'), (A. 5), ([^3)1 . (l3~T2]) hold, and 



that u and Uh are the solutions respectively of (jl.ip - (|1.2j) and (|3.1[) 

There are constants Kl,Kj > 0, 6 > such that if Ax £ (0,6) and At satisfies 
the CFL condition (|3?T2|) . th en in yu, 

—K(At 1 / w + Ax 1 / 5 ) < u — Uh < K(At 1 / A + Ax 1 / 2 ) for 7 e [0,1), 

-KiAt 1 / 10 + Ax 1 / 1 ") < u -U h < K(At 1 / 4 + Ax 1 / 4 ) for 7 e [1,2). 

Proof. Let us write the scheme (13. ip in abstract form (|4.1[) . To this end, set 
[u]t.x(s, y) = u(t + s,x + y) and divide p. lip by At to see that Q3.ip takes the 
form (|4.ip with 



(b n g' n Ja) ^K n Ja) 
S{h,t n ,X{3,r, [u\ tniX0 ) = sup < ■ r - ^ '— [u\ tri , Xfi (0,Xp - xp) 

rn,n—l/ \ \ 

fc^Jo) 1 

By its definition (|3.ip , monotonicity p. lip , and consistency (I3.9p . this scheme 
obviously satisfies assumptions (SI) - (S3) if the CFL condition (|3.12p holds. In 
particular, from (|33| and (pT5|) . (pO]) . ([57TD]> . wc find that 

ft, e) - E 2 (K, h, e) |^ (Afe _ 3 + Ax£ -i + (Aa;2 + ^-3), 7 £ ^ 2) . 
The result then follows from Theorem l4.2l and a minimization with respect to s. □ 

Remark 4.5. The error estimate is independent of 7 and robust in the sense that 
it applies to non-smooth solutions. 



5. New approximations of the non-local term 

In this section we derive direct approximations J£ [u] of the non-local integro term 
J a [u] appearing in (II. ip . As in pQ (cf. also [34]), the idea is to perform integration 
by parts to reduce the singularity of the measure. For the full discretization of (jl.ip 
along with convergence analysis, we refer to Section [3] 

We consider 3 cases separately: (i) ^ z ^ <1 v(dz) < 00, (ii) /i_ s i <1 \z\v(dz) < 00, 
and (iii) J^^j \z\ 2 v(dz) < 00. Note that in cases (i) and (ii) we can write the 
non-local operator in the form 

J a [<j)\{t,x) = I a [4>](t,x) -b a (x)Dcp, (5.1) 

where 

I a [4>]{t,x) ■= I U(t,x + V a (t,x,z))~Av{dz), 

J\z\>0 V ' 

b a (x) := [ if(t,x,z)v(dz), 

JQ<\z\<l 

for smooth bounded functions <fi. The reason is that I a \<j>\ and b a (x) are well-defined 
under assumptions (A. 2), (A. 3), (A. 4) if either (i) or (ii) holds. Furthermore, 
b a (x) will be bounded and x-Lipschitz. The term b a D(j) will be approximated by 
quadrature and upwind finite differences as in Appendix El leading to a first order 
method. Wc skip the standard details and focus on the non-local term /"[</>]. 

To simplify the presentation a bit, we will only consider the Cartesian x-grid 
{ x p}p — Ax7t N , but it is possible to consider unstructured non-degenerate families 
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of grids. On our grid we define a positive and 2nd order interpolation operator i^, 
i.e., an operator satisfying 

ih<fi(x) = wp{x)(j){xp) with wp{x) > 0, (5-2) 

/3ez N 

\E x [<l>](x)\ := \4>{x) - i h <t>{x)\ < K!&x 2 \D 2 cp\ , (5.3) 

for all x € M. N and where wp(x) > are basis functions satisfying w@{xp) = 5p p 
and Wp = 1. Linear and multi-linear interpolation satisfy these assumptions. 
Note that higher order interpolation is not monotone in general. 
We will also need the following monotone difference operators: 

S^ h (f>{r,y) = ±^{(j>{r ± Ax,y) - ^{r,y)} , (5.4) 

A rrik <f>(r,y) = — {<f>(r + k,y)-2^(r,y) + <f>(r-k,y)}, (5.5) 
for functions <p(r, y) on R x M. K for some K E N. For smooth </> we have 
\5± h <f> - d r <j>\ < i|0 rr |oAx, \A rr ^ - d 2 r cj>\ < l\dU\o\k\ 2 . 

5.1. Finite Levy measures. Assuming L^Hife) < oo, we approximate the 
term I a [(f>] defined in (HH]) by 

1% [4>] (t, x) = Q h [(i h (p)(t, x + r] a (t, x, z)) - <j>(t, x) 

where Qh denotes a positive quadrature rule on the z-grid {zp}/3 C M M with 
maximal grid spacing Az, satisfying 

Qh[4>]= X! ^^{zp) witn U P ^ °> 

l^oMI := IJX*M<fe) - QfcMI < K Q Az k *\D k Q<t>\ fv(dz), 

for smooth bounded functions </>, where Kq > and &q S N. Many quadrature 
methods satisfies these requirements, e.g., compound Newton-Cotes methods of 
order less than 9 and Gauss methods of arbitrary order. Note that the z-grid does 
not have to be a Cartesian grid. This method is at most 2nd order accurate because 

I a [0] = ™ + E^i; ■ + V a )]fv(dz) + E Q [<t>(; ■ + V a ) - 0], 

and it is monotone by construction, satisfying (|3.6[l and (|3.8p . The 0(Ax~ 1 ) term 
in (13.811 comes from the discretization of the b a term in (15. 



5.2. Unbounded Levy measures I. Now we assume that JL| <:1 |z|z^((fe) < oo, 
or more precisely that (A. 4') holds with 7 < 1. We consider the one-dimensional 
and multi-dimensional cases separately. 

5.2.1. One- dimensional case (M = 1). Now I a [4>] in (|5 . 1 [) takes the form 

I a [4>](t,x)= [ [(j)(t,x + n a {t 7 x,z)) - <t>{t,x))k{z)dz. 
We approximate this term by 

ih\m,x) 

oc 

H fcfcC^X*' x + rf(t, x, z n ))kf i n - 8~ h (i h <f>)(t, x + r] a (t, x, z- n ))k~ n 



71 = 



(5.6) 
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where z n = nAx, Sf h is defined in (|5.4p . the ^-interpolation ih satisfies (|5.2I) and 
(|5.3p . Moreover, 

\Kn--=i:Z + ^)dz {) -\CHC)dC, if Z>0. 

By (A.4') (M = 1 and 7 < 1), < f M k(z)dz < 00. 

To derive this approximation, the key idea is to perform integration by parts: 

I a [<t>](t,x)=(J )(^(t,x + ri a (t,x,z))-cP(t,x))k(z)dz 

d f° d 

— ((j>(t,x + rj a (t,x,z)))k(z)dz - / — (<j>{t, x + r] a (t, x, z)))k(z)dz, 



z J_ 00 dz 

for bounded C 1 functions <fi. Write I a [(j>] = I a ' + [tfi] + I a '~[4>], and use quadrature, 
finite differencing, and interpolation to proceed as follows: 



I a '+Mt,x) := / d z [<f>(t,x + r} a (t,x,z))]k(z)dz 
Jo 

00 

^^2d t [<f,{t,x + V a {t,x,z))]\ kl 



n=0 



E 4>(t, x + rj a (t, x, z n + Ax)) - 4>(t, x + r) a (t,x, z n )) L+ 
Ax h ' 



n=0 



E {ih<fi)(t, x + 77° ft, x, z n + Ax)) - (ih(f>){t, x + 77° ft, x, z n )) + 
Ax h < n ' 

n=0 

In a similar way we can discretize I a, ~[(j)] and (|5.6p follows. 
The approximation just proposed is consistent since 

I a [4>] (t,x)= 1% [4>] (t,x)+E Q + E FDM + Ej , 

where Eq, -Efdm, and Ei denote respectively the error contributions from the 
approximation of the integral (1st order), the difference approximation (up- winding, 
1st order), and the 2nd order interpolation. These terms can be estimated as 
follows: 



\E Q \ <Ax\d^{- + r ] a )\ / k(z)dz, 

Jm. 

\E FDM \ < ^Az^G + Olo [ k(z)dz, 
\E I \<2Ax\D 2 x cj)(-+T 1 a )\ { k[z)dz. 



The discretization (|5.6|) is also monotone satisfying ()3.6() and also (|3.8ll (when I£ 
replaces J£). To see this, note that ih4>(xp) = 0(cc^) and that by (A. 2) 77ft, x, 0) = 0. 
Hence we can reorganize the sum defining and write 

C + Mft,^) 

1 1 00 

— k+ <P(t,Xp) + — Y,( k h,n-1 - k h,n)( i h ( l ) )( t ^ + V a (t,Xp,Z n )) 



Ax 



n=l 



1 V 



Ax 

n—1 
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In a similar way 

1 x 

n—1 

Since k is increasing on (0, oo) and decreasing on (— oo,0), 

K h,n-1 ^ K h,n' 

and hence by (|5.2[) and J2f3 w P = 1> <|3.6|> and Ij3.8|) hold with 

i E »/»(^+^^.«))(<S-i-<w w )^ ' 



i£Z\{0} 

and k^ = ©(Ax 1 "^). The leading ©(A^ 1 ) term in (|3.8j) comes from discretizing 
the b a term in (|5.ip . 

5.2.2. Multi- dimensional case (M > 1). In this case we write I%[(/>] of (|5.ip in polar 
coordinates and propose the following approximation: 

£ <5+ ft [i h <f>(t, x + n a (t, x, r n y))} k h , n (y)dS y , (5.7) 

i^N 1 n=0 

where r n = nAi, cZSj, is the surface measure on the unit sphere in M M , 5± h is 
defined in (|5.4p . the ^-interpolation satisfies (|5.2p and (|5.3[) . Moreover, 

/•r n+ i ^ a ^oo 

kh,n(y) = / k(r,y)dr and k(r,y) = I k(sy)s M ~ 1 ds. 

J r n J r 

By assumption (A. 4') with 7 G (0, 1), < J °° k(r, y)dr < C < 00 for all |y| = 1. 

To derive this approximation we use polar coordinates and integrate by parts in 
the radial direction. Let be a bounded C 1 function, and set 

G a (t, x, z) := <f>(t, x + r] a (t, x, z)) - (j>(t, x). 

Then 

I a {(j)}{t,x) = [ G a (t,x,z)k(z)dz 

J& M \{0} 

poo 

[ / G a (t, x, ry)r M ~ 1 k(ry)dr] dS y 
|»|=i Jo 

d 

—G a (t, x, ry)k(r, y)dr] dS y , 



I 1/ 

J\y\ = l JO 



Or 

and (|5.7p follows by discretizing the inner integral as in Section lOTTl 
This is a consistent first order approximation of I a [4>] since 

I a [4>] (t,x) = 1% [4>] {t,x)+E Q + E FDM + E l , 

where Eq, -Efdm, Ej have the same meaning as in Section 15.2.11 and these terms 
can be estimated as follows: 



\E Q \ < Aa;|D^(. + O| M, 
I^fdmI < \&x \Dl<j,{- + T! a )\oM y 
\Ej\ < Ax\D 2 x <p(-+ V a )\ M k , 
where M k = Jj y|=1 J °° k(r, y)dr dS y . 
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The approximation I^[</>] is also monotone, and satisfies (|3.6p and (|3.8[) . This 
follows as in Section f5. 2. 1[ since I£[(f))(t, x^) can be written as 



1 

Ax 



X] [ k h,n-l{y) - k h ,n(y)] (ih<f>)(t,Xp +l] a (t,X0,r n y)) - (f>(t,Xp) 
l»l=ln=l 



where for fixed y, kh, n (y) is a decreasing function in n since k(r,y) decreasing in 
r. Moreover, Zo' n has a term like -^kh,o(y) = C(Ax~ 7 ) plus the leading 0(Ax^ 1 ) 
term which comes from the discretization of the b a term in (15.11) . 



5.3. Unbounded Levy measures II. We assume that j\ z \ Kl \z\ 2 v(z)dx < oo, or 
more precisely that (A. 4') hold with 7 £ [1, 2). In this case the decomposition (|5.ip 
is not valid. Again, we consider the one-dimensional and multi-dimensional cases 
separately. 

5.3.1. One- dimensional Levy process (M = 1). Now the nonlocal operator takes 
the form 

J a [<p] (t, x) = I [cj)(t,x + 'n a (t,x,z))~(t){t,x)-r] a (t,x,z)D(t)]k(z)dz, 

JR\{0} 

or, after two integrations by parts (more details are given below), 

J a [<j>](t,x) = J a ' + [(/>](t,x) + J a >-[<f>](t,x) -l a (t,x)D4>, (5.8) 
where b a (t,x) = d^T] a (t,x, z)k(z)dz, 

/■±oc 

J a ' ± [(/)}=± d 2 z [(/)(t,x + r] a (t,x,z))]k{z)dz, (5.9) 



k{z) 



J-oo /-oo k( y r ) dr dw ' for z < 
f f k(r)dr dw, for z > 0. 

By (A.T) (M = 1, 7 < 2), < £;(z) < C\z\ 1 ~"^e-^ A+E '>^ and jfe is integrable. 

Note that b a is bounded and rr-Lipschitz, and that b a D<\> can be discretized 
using quadratures and finite differences as in Appendix [A] This leads to a first 
order monotone (upwind) approximation - we skip the standard details. 

We propose the following approximation of J a ,:t [</>]: 

00 

J^M^x) = Y / ^Az[^(t,x + if(t,x,z n ))]~k^ n , (5.10) 

71=0 

where z n = nAx (not nAz!), A 22j az is defined in (|5.5p . the ^-interpolation ih 
satisfies (15.21) and (15.31). Moreover, 



fc+„ = / /c(z)dz and fc^„ = / fc(*)cte. 

The approximation (|5.10p can be derived from (|5.9p using quadrature, finite differ- 
encing, and interpolation. 

To obtain (|5.9p and (|5.8p . we integrate by parts twice: 



[<f>(t,x + T) a (t,x,z)) - (f>(t,x) - rf(t, x, z)D(f>]k(z)dz 

p OO 

[c/)(t,x + ri a (t,x,z))-<f>(t,x)-n a (t,x,z)D4>]( - / k(w)dw) 

J Z 

/oo poo 
d z [(f)(t,x + T] a (t,x,z))-T] a (t,x,z)D(f>]( J k{w)dw)dz 

d z [4>{t, x + r] a (t, x, z)) - n a (t, x, z)D<f>] ( - k) (z) 
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+ / d 2 z [(f>(t,x + r ] a (t,x 1 z)) -r ) a (t,x 1 z)D(f>]k(z)dz 
Jo 

/>oo poo 

+ 0+/ d 2 z [(j)(t,x + ri a (t,x 1 z))]k(z)dz- D<j> I d 2 rj a (t 7 x, z)k(z)dz. 



In view of this result and similar computations for the integral on (— oo, 0), (|5.8|) 
follows. These computations are rigorous if </3(t, x + 77), d z <j)(t, x + 77), d 2 (ft(t, x + ri) 
and 77, d z n, d 2 rj are z-integrable and bounded by e A ' 2 ' at infinity. 
The approximation is consistent and has the error expansion 

J*'* [0] (t, x) = J%>± [0] (t, x)+E± + E± DM + Ef, 

where Eq, £fdMj Ej have the same meaning as in Section 15.2.11 and these terms 
can be estimated as follows: 

\E%\ < Ax\d^{-+r] a )\o f k(z)dz, 

Jm 

I^fdmI < ^Az 2 \d 4 z ^- + v a )\o f Hz)dz, 

At 2 f 

\Ef\<A-^\Dld>(- + rj a )\ / k(z)dz. 

The proposed approximation is first order accurate if Az = Ax 1 ^ 2 , it is monotone 
satisfying (|3.6|) . and (|3.8p holds if Az = Ax 1 / 2 . These properties follow as in Section 
15.2.11 since J^ :± [4>] it, xg) can be written as 

^2 *M [(**0)(*> x + ^"(^ ^^tO) ^ ^)] 

1 ^ 

+ ~ 2 ^,n + Kn-l) [MX*. * + ^"C*. *±n)) ~ 0(*> ^)] > 

n— 1 

and, by convexity of k(z) on (0, cxd) and (— oo,0), 

Kn+l - 2l tn + Kn-l > for n > 1. 

Moreover, j^' n equals -£?(k£ Q — k'^ l + k^ Q — = 0(Ax 2 ^ 7 /Az 2 ) plus a 

C(Aa; _1 ) term from the discretization of the & a -term in (|5.8|) . When Az = Ax 1 / 2 
the leading term is the C(Aa; _1 ) term. 

5.3.2. Multi- dimensional Levy process (M > 1). Writing J a [(j)] in polar coordinates 
and performing two integrations by parts in the radial direction leads to 

J a [(j>}it,x) = J a [(/)}it,x) -b a it,x)D<P, (5.11) 
where b a it, x) = Li =1 J °° d 2 [r] a it, x, ry)] k(r, y)dr dS y and 

j a [<f>](t,x) = / ^ d^(t,x + r, a {t,x,ry))]k{r,y)drdS y , 

J\y\ = l JO 

oo />oo 

M-1 i 



k(s) = / r M - l kiry)drdw. 

J s J w 

By (A. 4') (7 < 2), k(r,y) < C'r 1_7 e _ ' A+e - ), ' and thus k is r-integrable uniformly 
in y. Note that b a is bounded and x-Lipschitz, and that b a D(j> can be discretized 
using quadrature and finite differencing as in Appendix |XJ This leads to a first 
order monotone (upwind) approximation - we skip the standard details. 
We propose the following approximation of J a [4>] : 

/oo 
x + 77 Q (i,a;,r„?/))] kh,n{y) dS Vl (5.12) 
,J/l =1 n=0 
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where r n = nAx (not nAz\), A rr ^ z is denned in (|5.5p . the x-interpolation ih 
satisfies (|5.2p and (|5.3p . and 

kh,n(y) = / k(r,y)dz. 



The approximation (|5.12[) follows from (|5.1ip by quadrature, finite differencing, 
and interpolation, and the derivation of (|5.11|) is rigorous provided the functions 
<fi(t, x+ij), D z tfi(t, x+rj), D 2 (f>(t, x+ij) and 77, D z -q, D 2 z r\ are z-integrable and bounded 
by e A ' z at infinity. 

The approximation is consistent and has the error expansion 

J a [4>] (t, x) = 3% [<f>] {t,x)+E Q + £ FDM + Ej , 

where Eq, _Efdm 7 Ei have the same meaning as in Section 15.3.11 and can be 
estimated as follows: 

\E Q \ < Ax\D 3 z <f>(- + r, a )\ M~ k , 
\E FDM \<^Az 2 \D^(- + V a )\ M h 

\ E i\<^\D 2 x 4>(-+V a )\oM~ k , 

where Mj, :— J\ y \—i k{r,y)dr dS y . Whenever Az = Ax 1 / 2 , this is a first order 
approximation. Moreover, the approximation is monotone satisfying <|3 - 6 j) and, 
whenever Az — Ax 1 ' 2 , it also satisfies (|3.8|) . This follows as in Section f5 . 2 . 1 1 since 
J^[<f>](t, Xa) can be written as an integral over {\y\ = 1} with integrand 

-J~2 [{ih<t>){t,Xp +?i a {t,xp,r_ 1 y)) - (j)(t, x fj )}~k hfi {y) 



Az 2 ^ 

n=l 



^ (k h}n+1 (y) - 2kh, n (y) + kh,n-i(y)) 



\{i>h4>){t,- 



,Xp + V a (t,Xp,r n y)) - (f)(t,Xp)] 
Furthermore, for each fixed y, fc(r, y) is convex on (0, 00) and thus 

kh,n+\{y) - 2kh, n (y) + kh,n-i(y) > for n > 1 and \y\ = 1. 
Remark 5.1. 

a. (Order of schemes) In general, our discretizations of the non-local term in 
are at most first order accurate. In the case 7 e [1,2), a first order rate is obtained 
by choosing Az — Ax 1 / 2 . Higher order discretizations can be derived using higher 
order quadrature and interpolation rules, but the resulting discretizations are not 
monotone in general. On the other hand, if 77 = z, then interpolation is not needed 
and consequently the monotone discretizations of Section ^. 3l arc 2nd order accurate. 



b. (Remaining discretizations) To obtain fully discrete schemes it remains to 
discretize the various terms involving b a D(jj, for example by quadrature and finite 
differencing, cf. Appendix [X] In applications, the densities k and k can often be 
explicitly calculated, e.g., using incomplete gamma functions as in 1 . Otherwise 
these quantities also have to be computed by quadrature. Furthermore, regarding 
Sections 15.2.21 and 15.3.21 it also remains to discretize the surface integral in y. This 
discretization does not pose any problems, neither numerically nor in the analysis, 
as long as positive quadratures are used. The details are left to the reader. 
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c. (Increasing efficiency) From a practical point of view in terms computational 
efficiency, quadratures should be implemented using FFT. This is standard and we 
refer to, e.g., [22] for the details. 

d. (Generalization I) The above approximations (with obvious modifications) 
also apply to integral terms of the type 

J a [<l)}{t,x) = V / (<j>(t,x + Ttf(t,x,z))-</>(t,x) -r]f{t,x,z)D<j)(t,x))ki{z)dz. 
i=1 MM 

Such terms appear in M-dimensional Levy models based on M independent Poisson 
random measures coming from one-dimensional Levy processes. This is a rich class 
of models with many applications. For more information and analysis of such 
models we refer to the book [38 . 

e. (Generalization II) With obvious modifications, our approximations also apply 
to linear and non-linear equations involving the fractional Laplace operator 

f u(x + z) - u(x) - zDu(x) 
(~A) a u{x) = c a -i v ^ K -^dz, 6(0,1), 

J\z\>0 \ z \ 

where x, z G and c a is a constant, in which case the Levy measure takes the 
formz/(dz) = \z\- N - 2a dz. This measure satisfies |(A.3)| except for the "exponential 
decay at infinity" requirement. It is straightforward to recast the entire theory to 
allow for a fractional Laplace setting where assumption (A. 3) is replaced by 

\z\ 2 A 1 v(dz) < oo. 

z|>0 

6. Error estimates for a switching system approximation 

In this section we obtain error estimates for a switching system approximation 
of (|1.1|) - (|1.2|) . This result, which has independent interest, plays a crucial role in 
the proof of Theorem 14.21 in Section [7J 

The switching system will be written as 

F i (t,x,v,d t v i ,Dv i ,D 2 v i ,v, i (t r ))=0 in Q T , i G {1,2, ,m}, (6.1) 

v(0,x) = (g(x),...,g(x)) in R N , (6.2) 
where v = (v\, . . . , v m ) is in W 1 and for sets Ai such that UiAi = A, 
Fi(t,x,r,p t ,p x ,X,<p(-)) 

= max{p t + sup [£ a (t,x,ri,p x ,X) - J a [4>](t,x)];ri- M*r}, 

aeAi 

C a {t,x,r,p,X) := -tr(a a {t,x)X) - b a (t,x) ■ p + c a {t,x)r - f a (t,x), 
M l r = minjVj + k}, k > 0, 

for (t,x,r,p t ,p x ,X) e M x R N x W m x R x R N x § w and any smooth bounded 
function 4>. The operator J Q [</>] is defined below (|l.l|l . 

In the pure PDE case, such approximations have been studied in, e.g., [I7 | [23] , [8] . 
Here we extend the error estimates of [8] to non-local Bellman equations. In a 
complimentary article |12j , we develop a viscosity solution theory covering switching 
systems like (|6.ip . We refer to that paper for the precise definition of viscosity 



solutions and proofs of the associated results utilized herein. If assumptions (A.l) 



(A. 4) of Section [3] hold, then we have the following well posedness result |12j 



Proposition 6.1. Assume that conditions (A.l)- (A. 4) hold. There exists a unique 
viscosity solution v of (|6.1[) - (|6.2[) , satisfying \v\± < C for a constant C depending 
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only on T and K from (A.l) - {A.3)\ Furthermore, if w\,W2 are respectively vis- 
cosity sub and supers olutions of (|6.1|) satisfying wi(0, •) < ^(O, •), then w\ < w%. 

Before we continue, we need the following remark. 

Remark 6.1. The functions a a ,b a ,c a , f a , r] a are only defined for times t G [0, T]. 
But they can be easily extended to times [— r, T + r] for any r > in such a way 



that (A.l) - (A. 3) still hold. In view of Proposition 16.11 we can then solve the 
initial value problem up to time T + r or, by using a translation in time, we may 
start from time — r. We will use these facts several times below. 

By equi-continuity and the Arzela-Ascoli theorem it easily follows that the each 
component of the solution of l|6.ip - (|6.2p converges locally uniformly to the solution 
of p.ip - (|1.2p as k — > 0. To derive an error estimate we use Krylov's method of 
shaking the coefficients coupled with an idea of P.-L. Lions as in [5]. We need the 
following auxiliary system 

Ff(t,x,v\d t v!,DvlD 2 vlvt(t,-)) = in Q T+£ 2, ie{l,...,m}, (6.3) 

v^0,x) = (g(x),...,g(x)) in R N , 

where v e = (v\, . . . , v^) and 
F-(t,x,r 1 p t ,p x ,X, 4>{-)) = max {p t + sup (C a (t + s,x + e,n,p x , X) 

a€Ai;\e\<e;0<s<e 2 

- J a {t + s,x + e)(j));n - M l r) . 

The operators £, J, and M are as previously defined. 

Note that we have used the extension of the data mentioned in Remark 16. II By 
regularity and continuous dependence results from |12j we have 



Proposition 6.2. Assume that (A.l)- (A. 4) hold. There exists a unique viscosity 



solution v e : Qt+c 2 — > K of (|6.3p satisfying 

Mi + V-f|o <C, 

e 

where C depends T and K . Furthermore, if w\ and u>2 are respectively sub and 
supers olutions of (|6.3p satisfying Wi(0, •) < 102(0, ■), then wi < W2- 

We are now in a position to prove the following main result of this section: 



Theorem 6.3. Assume that (A.l)- (A. 4) hold. If u andv are respectively viscosity 
solutions of (|l.ip - (|1.2p and (|6.ip - (|6.2p . then for sufficiently small k, 

0<Vi — u<Cki, ie{l,...,m}, 
where C depends only on K and T. 

Proof. Since w = (u, . . . , u) is a viscosity subsolution of (|6.ip . the first inequality 
u < Vi follows from the comparison principle. 

The second inequality will be obtained in the following. Since v e is the viscosity 
solution of (|6.3p . it follows that 

d t vf+ sup (C a (t + s,x + e,v £ l (t,x),Dv e l ,D 2 v e l )-J a (t + s,x + e)vt) <0 

aeAi 

in Qt+<l 2 in the viscosity sense, i — 1, . . . , to. After a change of variable, we conclude 
that for every < s < e 2 and \e\ < e, v e (t — s,x — e) is a viscosity subsolution of 
the uncoupled system 

d t wl + sup (C a (t, x, w\,Dwl, D 2 wt) - J a {t, x)w-) = in Q e T , (6.4) 

aeAi 

where :— (e 2 , T) x K . Now set v e :— v e * p e , where p e is the mollifier defined 
in p.3p . A Riemann-sum approximation shows that this function is the limit of 
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convex combinations of viscosity subsolutions v(t — s,x — e) of the convex system 
(|6.4p . Hence v t is also a viscosity subsolution of (|6.4|) (see the appendix of [25] for 
more details). On the other hand, since i> £ is a continuous subsolution of (|6.3|) . 

v\ < min Vj + k in Q T . e 2, i G {1, . . . , m). 

It follows that maxi — mintij(i,x) < k in Qt+c 2 ; an d therefore 

k e -wj|o<fc, i,i e {l,...,m}. 
Then, by the definition and properties of v e , 

k k 
\d t v tt -d t v ej \o<C— and \D n v tt - D n v ej | < C— , 

for n e N, i,j 6 {1,..., m}, where C depends only on p, T, and K. For e < 1, it 
follows that 

\d t v ej + sup (C a (t, (t,x),Dv ej ,D 2 v ej ) - J a (t,x)v tj ) 

aeAi 

k 

-d t v ti - sup (C a {t, x,v t i(t,x),Dvti,D 2 v <ii ) - J a (t, x)v ti )\ < C— 
and, since v e is subsolution of (16.41). 



k 

d t v ti + sup x,v ei (t,x),Dv ti ,D 2 v ei ) - J a (t,x)v ti ) < C— in Q e T , 

aeA e 

where the constant C depends on p, T, and K . From this inequality it is easy to 
see that v e i — te Kt C^ is a subsolution of (jl.ip restricted to Qt- Hence, by the 
comparison principle, 

v ei ~ u < e Kt (\v el (e 2 , •) - u(e 2 , -)|o + Ct^j in Q e T , i € {1, . . . , m}. 

By regularity in time, •) — Vi(t, -)|o < (|w|i + |fi|i)e, and by Proposition 16.21 
and properties of mollification we conclude that 

k 

Vi - u < Vi - v ti +v ei -u< C(e + -=■) in Q T , i G {1, . . . , m}. 
Now the theorem follows by minimizing with respect to e. □ 
7. The Proof of Theorem 14.21 



To prove Theorem 14.21 we will use different arguments for the upper and lower 
bounds. The upper bound, part (a), is the "easy" part, and it is essentially a 
reformulation of the general upper bound established in [53]. We skip the details, 
and prove only part (b) which is a new result. 



Proof of Theorem \4-%\ (b). Without loss of generality we will assume that A is finite: 

A = {ax, a 2 , ■ ■ ■ , a m }. 

The proof of this statement is similar to the one given in [S] in the pure PDE case 
and relies on assumption (A.l) Now we follow [H] and use a switching system 
approximation to construct approximate supersolutions of (jl.ip which are point- 
wise minima of smooth functions and approximates the viscosity solution of (11. ip - 
([L2]) , Consider 

Ft(t,x,v\d t vlDvi,D 2 viy i {t,-)) =0 inQ T+2e2 (7.1) 

v e (0,x) = v Q (x) in R N , 
where v e = (uf, . . .,v e m ), v Q = (g, . . .,g), and 
Ff(t,x,r,p t ,p x ,X, cf)(t, •)) 
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= max {p t + min (C ai (t + s — e 2 , x + e, r, p t , p x . X) 

1 0<s<e 2 ,|e|<e V V ' 

- J a ' (t + s - e 2 , x + e)4>) ; r< - M l r) , 

where £, J and M. are denned below (|6.ip in Section [5] This new problem is well- 
posed and each component of the solution of this switching system will converge to 
the viscosity solution of (jl.ip as k, e — » 0: 



Lemma 7.1. Assume that conditions (A.l) - (A. 4-) hold. There exists a unique 
solution v e of (|7.1[) satisfying 

\v e \\ < K, max |w| — v e ; \q < k, and for k small, max \u — vf\ < C(e + k^), 

i,3 i 



where K,C only depend on T and K from (A. 2) - (A. 4-) 



Proof. From [12] we have the existence and uniqueness of a viscosity solution, and 
moreover the uniform bounds 

\v e \i<K and |u e -v°\ < Ce, 

where v° is the unique viscosity solution of (|7.ip corresponding to e — 0. The 



last inequality in the lemma now follows since \u — u?|o < Ck~s by Theorem | 
To second inequality follows since arguing as in the proof of Theorem 16.31 leads to 
< maxi v| — minj Vj < k in Qt+2<l 2 ■ D 

Next we time-shift and mollify v e . For i = 1, . . . , m, set 

vl(t,x) := vl(t + e 2 , x), v ei (t,x) := p t *vl(t,x), 

where p is defined in (|1.3|) . Note that supp(p £ ) C (0,e 2 ) x B(0,s) and that the 
functions v e , v e , v e are well-defined respectively on Qt+2s 2 > (-£ 2 ,T + e 2 ] x M. N , 
Q T+e 2. By Lemma 1 7. II and properties of mollifiers, 

|^ e |i<^, \v'-v e \ a <Ke, 

max \v ei - v e j\ < C(k + e) in Q T +e 2 , (7-2) 
max |u — ti £j i| < C(e + /c3 ) in Qt, 



where C depends only on p and K, T from (A. 2) (A. 4) A supersolution of 
can now be produced by setting 



w 



Lemma 7.2. Assume that conditions (A.l) — (A. 4) hold and e < (8 supjufji) 1 k 



For every (t,x) G Qt, if j '■= argmin^g^t, a;), 

d t v ej + C a > (t 1 x 1 v e] (t,x),Dv ej {t 1 x) 1 D 2 v ej {t,x)) ~ J a i{t,x)v ej > 0. (7.3) 

We postpone the proof of this lemma. From this lemma it follows that w is an 
approximate supersolution to the scheme (|4.ip when e < (8 supjti|]i) _1 fc: 

S(h,t,x,w(t,x),[w] t , x ) > -E 2 (K,h,e) in Q+ , (7.4) 

where K conies from Lemma 17.11 To see this, let (t, x) € Qt and set j :— 
argmin^ci^, x). At (<, x), w(t,x) — v t j(t,x) and w < v e j in Qh- Hence (SI) 
implies that 

S(h,t,x,w(t,x), [w]t :X ) > S(h,t,x,v ej (t,x), [v ej ] t . x ). 
By consistency (S3) (ii) we have 
S(h,t,x,v 6j (t,x), [V e j]t,x) 

> d t v ej + F(t,x,v ej (t,x),Dv e j,D 2 v ej ,v e j(t, •)) - E 2 (K,h,e), 
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> dtVy+£ a > (t,x,v ej (t,x),Dv ej (t,x),D 2 v ej (t,x)) ~ J a > (t, x)v ej , - E 2 (K, h, e), 

and (|7.4[) then follows from Lemma [731 

To derive the lower bound on the error Uh ~ u, we take e = (8 supjiifji)" 1 ^ and 
use (|7.4p and comparison Lemma 14.11 to get 

u h - w < e^\(g h - w(0, ■))+{ + 2te^E 2 (K, h, e) in Q h . 

By ([7T2]) . \w-u\< C(e + k + ki), and hence 

u h -u<e^\(g h -w(0,-))+\ + 2te^E 2 (K,h,e)+C(e + k + k?) in Q h , 

possibly with a new constant C. Since e = Ck, the proof is complete by minimizing 
the right hand side with respect to e. □ 

Proof of Lemma \ 7. 2\ We begin by fixing {t, x) £ Qt and set j = argmin^ v e i(t, x). 
Then 

v e j(t, x) — A4 J v e (t, x) = max {v e j(t, x) — v e i — k\ < —k. 

Therefore, by the Holder continuity of v e and basic properties of mollifiers, 

v € Jt,x) - M j v £ (t,x) < -k + 2max[^ e ]i2e 

J i 

and 

v e Js,y) - M 3 v e (s,y) < -k + 2max[<]i(2e + \x - y\ + \t - s|i), 

for all (s,y) G Qt- Consequently, if \x — y\ < e, \t— s\ < e 2 , and e < (8supjw|]i) _1 /c, 
then 

v){s,y)-M'v £ { Sl y) < 0. (7.5) 

To continue we need the following remark. Let u 1 ,. . . ,u k be functions satisfy- 
ing (|7.5p at (t,x), then any linear combination u A = 53j=i ^i u% with A; > and 
= 1 , also satisfies (|7.5j) at (t,x). If, in addition, u 1 ,. . . ,u k are supersolu- 

tions of 

max {d t Uj + C a ' (i, x, u ] ,Du j ,D 2 u J ) - J a ' [uj](t, x);u 3 - M j u) = 0, (7.6) 

then in view of (|7.5[) they are also supersolutions of the linear equation 

d t Uj + £ a > (t, x, u^Du^D 2 ^) - J a i [ Uj ]{t 7 x)=0 (7.7) 

at (t, x). An easy adaptation of the proof of Lemma 6.3 in [29 then shows that u x 
is also a viscosity supersolution of (|7.7[) at (t, x) . 

A change of variables reveals that {v £ (- — s, • — e)}^ , < s < e 2 , |e| < e, is 

a family of supersolutions of (|7.6| with [v e — M^v e ](t — s,x — e) < 0. We note 
that by approximating the function v e j by a Riemann sum, we see that it is the 
limit of convex combinations of v e (- — s, • — e). In view of the above remark, these 
convex combinations are supersolutions of (|7.6p . and hence by the stability result 
for viscosity supersolutions, so is the limit v e j. Finally, since this function is smooth 
it is also a classical supersolution of (17. 6|) and hence (|7.3j) holds. □ 
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Appendix A. An example of a monotone discretization of L a 

Let {ei}fL 1 be the standard basis of R N and the ij-th element of the matrix 
a. Kushner and Dupuis [33] suggest the following discretization of L a in (jl.ljl : 

N 

i=l ijtj 

where b + — max{6, 0}, b~ = (— b) + , and 
8f4>(x)=±-±-{cl)(x±e i Ax)-cl>(x)}, 

Au(f)(x) = '^^{ < l ) ( x + eiAx) - 2<j>{x) + <j>(x - e^Ax)}, 

A tj<t>{ x ) = ^-^i 2 ^) + <f>( x + e * Ax + e j Ax ) + <f>( x ~ e * Ax ~ e j Ax )} 

— ^ 2 {4>(x + eiAx) + 4>(x — eiAx) + <fi(x + ejAx) + <fi(x — e^-Aa;)} 

A^j(f)(x) — 2/± x 2 {^( x ) + + ^Ax - ejAx) + <f>(x - eiAx + e^Aa;)} 

— ^ 2 {4>(x + eiAx) + <fi(x — eiAx) + <fi(x + ejAx) + <fi(x — ejAx)}. 

By Taylor expansion it is easy to check that the truncation error is given by (|3.3p 
Moreover Lh can be written in the form (|3.5p with 

- 1 —ekmi 1 //; (/ " rl 



h,0,/3±ei 



Aa 



ja.n _ a jj fti x ) ja.n _ a jj fti x ) ■ , ■ 

i h,fi,l3+e i ±e j - ^ x 2 ' l h,l3,f)-e i ±e j ~ ^2 ' 1 T 3, 

and Ifr'pp — otherwise. This approximation is monotone if l^'^ ^ > for all 
a, (3,(3,71 and Aa; > 0, which happens to be the case, e.g., if a is diagonally domi- 
nant: 

af { (t,x) — 2_. \ a ?j{t? x )\ > in Qt, for each a e A. 
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